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Abstract. The ionizing spectral energy distribution of quasars is known 
to exhibit a steepening of the distribution short-ward of 1050A. It has 
been assumed that this change of power-law index from —1 to —2 is in- 
trinsic to the quasars. We study an alternative interpretation, in which 
a tenuous absorption screen is responsible for the change of index. A 
successful spectral fit is obtained using a total hydrogen density distribu- 
tion decreasing with redshift as ~ (1 + zp with £ of the order —3.5 if we 
suppose that the ionization state of the intervening gas is controlled by 
the line-of-sight quasar. This interpretation is not compatible, however, 
with the relative weakness of the observed quasar intensities with respect 
to that of the local background metagalactic radiation. Therefore, our 
interpretation in its present form cannot be sustained. 



1. Introduction 

The ionizing spectral energy distribution (hereafter ISED) of nearby active galac- 
tic nuclei cannot be observed directly due to the galactic absorption beyond the 
Lyman limit. Owing to the redshift effect, however, we can get a glimpse of 
the ISED from the spectra of very distant quasars. Most studies on the subject 
indicate a steepening of the distribution at wavelengths shorter than 1000 A 
(see O'Brien et al. 1988) and references therein). The most extensive work on 
the subject has been performed by Zheng et al. (1997, ZKTGD) using HST 
archived data. They found that the power-law index (F u oc u a ) steepens from 
~ —1 for A > 1050 A to ~ —2 for shorter wavelengths. 

Korista, Ferland & Baldwin (1997) pointed out that such a steep slope for 
the ISED would imply an insufficient number of photons beyond hv > 54.4 eV 
to account for the observed luminosities of the high excitation emission lines. 

We report preliminary results of our attempt to find an alternative expla- 
nation for the break, one that is extrinsic to the quasar ISEDs. We postulate 
the existence of a very tenuous absorption gas component pervading the local 
universe and proceed to study the characteristics it must possess to reproduce 
the observed break, and its place and role within our IGM views. 
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2. Data simulation to a first order 

The composite spectrum of ZKTGD was constructed by merging 284 spectra of 
101 quasars taken with FOS of the Hubble telescope. Although the wavelength 
coverage of individual spectra was relatively narrow, the distribution in quasar 
redshifts allowed ZKTGD to cover the wavelength range of 310-3000 A in the 
quasar rest-frame. Before merging their spectra, each spectrum was corrected 
for the Lyman valley and Lya forest absorption by calculating the appropriate 
transmission curve using the scheme developed by M0ller & Jakobsen (1990). 

Our aim is to reproduce the change of ISED index by introducing a pre- 
viously unknown intergalactic absorption component. This component must 
consist of clumps which are sufficiently optically thin (and numerous) that they 
do not show up as individually detectable Lya absorption lines in the archived 
FOS spectra, hence N H o <C 10 12 cm~ 2 . In this preliminary work, we did not 
fully simulate the observed data on a pixel by pixel basis, rather we defined a 
single spectrograph wavelength window (3000 A to 1300 A) which we considered 
typical of the wavelength coverage of individual quasars and then redshifted this 
window in locked steps. We defined a source array, Fj, to represent the intrinsic 
SED of all quasars. This SED consisted of a power-law of index a. All arrays 
were evenly binned in LogX and the wavelength of each array element j was 
such that LogXj = LogXuAX — j x u, with u = Log(l + 5z Q ) and Xmax being 
the highest rest- wavelength considered. The advantage of this scheme is that 
the source array Fj can be redshifted by simple translation of the pointers j 
provided the redshift is a multiple of 5z Q . This is not a serious limitation if one 
chooses a sufficiently small Sz Q (we used 0.01). The procedure followed consisted 
in translating the source array by an integer number of bins k and then add it 
to our averaging buffer: Bj = Bj + WkFjTj(zn, Xj) where Tj is the transmis- 
sion function evaluated at Xj for the k th array spectrum (corresponding to a 
quasar's redshift Zq = k5z Q ), and Wk is the global weight given to the k th array 
before summing. These 'uniform' weights were based on the dispersion vs. rest- 
frame wavelength relation displayed in the Fig. 6 of ZKTGD. More specifically, 
w k = s k 2 where Sk = A^ ' 71 with Afc c = 2150(1 + Zq) A, the central wavelength 
of redshifted array k. The lowest and highest redshifts considered were 0.33 
and 3.5, respectively, which are the values quoted by ZKTGD for their data set. 
We uniformly covered this redshift interval by considering all integer values of 
k between these two redshift limits. 

3. The transmission function 

For each quasar rest-frame wavelength bin j, we calculated the transmitted 
intensity = Ix^Tx^ = I\ j e~ T ^ by integrating the opacity along the line-of- 
sight to a quasar of redshift zq 

r(Xj) = E!£o Jo Q °i{ A) "M')£ dz (1) 

where Xj is the quasar rest-frame wavelength for bin j and njjo(z) the inter- 
galactic neutral hydrogen density. The summation was carried out over the 
different opacity sources: photoionization (i = 0) and line absorption from the 
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Lyman series of hydrogen (1 < i < 10). Although our code could include up to 
40 levels, we found that considering only the 10 lowest proved to be adequate. 
We adopted a fiducial velocity dispersion b of 30 km/s and assumed a simple 
line Gaussian profile for Cj of the lines. Apart from the uniform density absorp- 
tion case, we also implemented the treatment of a clumpy medium (described 
later) which was used to ensure that we could reproduce the transmission curve 
of M0ller h, Jakobsen (1990) or ZKTGD using parameters appropriate to the 
Lyo; forest. The self imposed constraint that A^o <C 10 12 cm" 2 for the clumps 
ensures that we remain in the linear regime of the curve of growth. Therefore, 
using a uniformly distributed gas or a clumpy medium will be equivalent. The 
use initially of an homogeneously distributed gas will save us from having to 
characterize the dumpiness using poorly defined parameters. 

For the dependence of the total gas density on redshift, let us consider a 
simple power-law tih oc (1 + z) s . What is lacking is the neutral fraction of nn 
which we will assume to be small and set by equilibrium between photoionization 
and recombination. In that case, we expect n H o to scale proportionally to 
the (total) density square, and inversely proportionally to the ambient ionizing 
radiation. We will consider negligible the cumulative opacity of the absorber (as 
confirmed a posteriori). With those considerations in mind, we obtain that n H o 
should scale as follows with z: 

n H0 {z)=n% ) (l + z) 2 ^T- 1 (z){l+Lo)- 1 (2) 

where the upperscript throughout the text denotes quantities evaluated at 
z = 0. In the above expression, we take into account the increase of density 
with the cosmological scale factor: n^/n^ oc (1 + z) 3 . The photoionization rate 
coefficient T due to the metagalactic ionizing radiation was parameterized as 
follows 

r(z) = (1 + z )0.73 e -0.526((*-2.3) 2 -2.3 2 )) (3) 

This is a renormalized version of the parametric fit of Haardt & Madau (1996) 
who modeled the metagalactic ionizing radiation (J^ CK ) taking into account the 
population evolution of quasars and their luminosity with redshift as well as self- 
emission of the clumpy intergalactic ionized medium. Note that T(z = 0) = 1 
because of our renormalization. 

The factor (1 + u;) _1 in equation 2 is the same as that defined by Bajtlik, 
Duncan & Ostriker (1988) and represents a decrease in the neutral fraction of 
hydrogen towards the line-of-sight quasar as a result of the ionizing radiation 
from the quasar. In this expression, the "proximity" term oj represents the 
geometrical dilution of the line-of-sight quasar ionizing flux with respect to the 
background ionizing radiation: 

- = m 2 < 4 > 

where ri is the luminosity distance of the quasar from a parcel of gas at z. 
if one of our free parameters and correspond to the distance from the quasar 
where its mean intensity at uq equals that of the background radiation J„ ■ 
At such distance, the neutral fraction x = n H o/riH is halved due to the quasar 
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flux. r 1 / 2 can be expressed as follows 

„o 

1/2 

with 



ri/2 = r? /2 J(l + ^)«- 1 r-i(*) (5) 




(6) 

where a is the power-law index of the quasar spectral luminosity distribution 
(L® oc v a ) and uq is the hydrogen photoionization threshold frequency. As 
ISED index, we adopt a = —0.86 which is the average index for the radio-quiet 
subsample of ZKTGD (above 1050A). In our data simulation, to describe the 
behavior of n H o with redshift, we employ only two free parameters: £ and r^ 2 . 

In equation 5, the quasar redshift as seen from the absorbing gas at z is z'q = 

(1 + zq)/(1 + z) - 1 while the Hubble constant is H' = H ^tt A + n M (l + z)' 6 
for a flat universe. We assumed the popular ACDM universe with $7a = 0.7 
and Q.m = 0.3. To compute the quasar luminosity distance 7\£,(zq), we used the 
parametric fit of Pen (1999). The light travel distance used to derive dl/dz was 
obtained by integration of the Freedmann equation. The adopted value for the 
Hubble constant is Hq = 67km/s/Mpc. 



4. Results 

Producing a break near 1050 A followed by a steepened power-law of index ~ 
—2 is fairly easy to achieve if one disregards any evolution of the background 
radiation (achieved by setting T = 1 in Eqn. 2) and adopt the parameters 
£ = -3.5, r° /2 = 3200 Mpc and n° Ho = 1.9 x 10 _11 cm _2 . This zeroth order 

model is illustrated by the short-dashed line in Fig. 1 (Model I) and is found to 
be relatively successful in producing a flat soft ISED. We should point out that 
photoionization by the line-of-sight quasar, as represented by the parameter 
r°/ 2 , is an essential ingredient to any successful fit. Actually, if we remove 
the uj term, we consistently obtain a sharp discontinuity at Lya, as illustrated 
by the short-dashed dotted line of Model II in Fig. 1. Including line-of-sight 
quasar photoionization smooths out completely the discontinuity provided r^ 2 
is sufficiently large. 

Model I presents, however, the following deficiency. We have compelling 
evidence that the background radiation increases strongly with z (c.f. Shull et 
al. 1999; Dave & Tripp 2001) and references therein. For this reason, varying 
r(z) in equation 2 using the function proposed by Haardt & Madau (1996) is 
preferable and certainly justified. Models III-VI are such models which span 
the range —5 to —2 in £. In these models, the parameters n° Ho and ry 2 were 
varied in such a way that the different curves superposed as much as possible 
one another within the interval 500-1000 A. All of our model parameters are 
listed in Table 1. We consider Model IV (solid line) to be the one providing the 
best overall fit. Notice that a broad hump near 500 A is apparent in all four 
models which consider J^ K evolution. The hump is only present in the case 
of the ISED derived from the radio-quiet subsample, which is the one plotted in 
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Figure 1. The composite spectrum of radio-quiet quasars con- 
structed by Zheng et al. (1997) (dotted line). The straight short- 
dashed line represents a power-law fit (a = —0.86) above 1050A to the 
composite continuum underlying the emission lines. The solid line is 
our favored Model IV (£ = —4, see Table 1). Three other models using 
different values of £ are represented by the long-dashed lines. From the 
top to the bottom curve, £ is decreasing from —2 to —5. The horizontal 
short-dashed line is Model I which does not include redshift evolution 
of the background radiation. Removing from Model I photoionization 
by line-of-sight quasar leads to Model II which is represented by the 
short-dashed dotted line. 
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Figure 2. The transmission curve as a function of observed wave- 
length for our best fit Model IV in the case of the quasar 0308—193 at 
z = 2.756 (solid line). Outside the FOS spectral window (1315-2400A), 
the dotted line shows the continuation of the calculated transmission 
curve. Starting from the right, the first, second and third descending 
ramps are contributed in large part by the opacity of Lya, Ly/3 and 
Ly7 lines, respectively. The small contribution of photoionization to 
the opacity is shown by the dotted line labeled phot, absorption. The 
short-dashed line represents the transmission curve for the case of a 
lower redshift quasar at z = 0.33. 
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Table 1. Model parameters 



Model 
(ID) 


n H0 
(cm -3 ) 


r l/2 

(Mpc) 




T evol. 
(Yes/No) 


I 


1.9 x ltr 11 


3200 


-3.5 


N 


II 


1.9 x 1CT 11 





-3.5 


N 


III 


7.4 x 10 -11 


3200 


-5.0 


Y 


IV 


4.3 x 1CT 11 


3200 


-4.0 


Y 


V 


3.2 x KT 11 


3000 


-3.0 


Y 


VI 


2.4 x KT 11 


2900 


-2.0 


Y 



Fig. 1. This hump is more apparent than in the ZKTGD paper because of our 
slightly lower wavelength limit at the UV end. To ascertain the reality of the 
hump would require a more extensive sample than the current. 



5. discussion 

In appearance successful, our alternative explanation of the steepening of the 
mean quasar ISED cannot be sustained for the following reasons: 

• A In our calculations, we have treated r®, 2 as a free parameter. There- 
fore, although we adopted a valid prescription for the variation of the 
background with z (eq. 3), we have not considered the issue of the absolute 
intensity of the background with respect to the observed quasar intensi- 
ties. It is straightforward to estimate u = (F u /4ttJ^ ck ) 2 at z = if we 

adopt an observed quasar flux of F^'° ~ 2 x 10~ 14 ergcm _2 s _1 near 1000 A, 
which is typical of the quasar sample of ZKTGD. The ionization rate coef- 
ficient of Haardt & Madau (1996) of 4 x 10~ 14 s^ 1 at z = translates into a 
mean intensity of J^ CK, ° ~ 1.3 x~ 23 ergs _1 sterad _1 cm _2 Hz _1 , hence u> ~ 
5.5 x 10 _27 /1.6 x 10~ 22 = 3.3 x 10~ 5 (using F VQ = 912 2 x 3.3 x 10- 19 F A() ). 
This contradicts our requirement that the proximity effect must extends 
up to z = for all quasars whose distance is less than 3 Gpc (the approx- 
imate value of Ty 2 ). This estimate of the local value of uj clearly shows 
that the background mean intensity at low z is much stronger than that of 
the quasar flux reaching us. Another way of illustrating this is to compare 
our values of r 1 / 2 with the high redshift value empirically determined of 
r i/2 ~ 8 Mpc by Bajtlik et al. (1988). If we translate this estimate to 
z = by taking into account the variation of J„ (z), we obtain that the 
Bajtlik et al. local value is r° , 2 w \/40 x 8 = 50 Mpc which is sixty times 
smaller than the value of 3200 Mpc proposed in Table 1. This discrepancy 
in the values of r?, 2 (which relates the quasar intensity to the relative 

strength of J^ CK, °) effectively rules out our simple model. 
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• B — As apparent in Fig. 1, our synthetic ISED rolls down steeply at the 
lower wavelength limit. If we calculate a similar transmitted ISED but 
for a reduced range in quasar redshifts (reducing the highest redshift to 
1, say), the roll over would be shifted to higher wavelength, a feature not 
seen in any individual quasar spectrum. 

• C The self-imposed requirement that individual absorbers or clumps be 
of columns Njjo <C 10 12 cm" 2 implies the existence of quite a large number 
of clumps along the line-of-sight. For a column of 10 10 cm -2 , about 6 x 10 7 
clumps per unit redshift would be required with an approximate individual 
mass of about one solar masses (assuming a neutral hydrogen fraction of 
1CT 3 and a clump hydrogen density of order 2 x lCT 4 /cc). 

In the case of point C, the number of clumps and their masses is not an 
unsurmountable objection, however, given the level of speculation characterizing 
current ideas about the nature of the dark matter. Point B above can also be 
resolved by having the proposed absorption component disappear at current 
Epoch (as a result of Compton heating to T > 10 6 K, say). The strongest 
objection is given in which definitely rules out the interplay between quasar 
and IGM photoionization as the correct physical mechanism for causing the 
smoothing out of the expected sharp break near 1000A (as the one seen in Model 
II of Fig. 1). The next step might consist in looking for a different distribution 
of neutral H than that given in Eqn. 2 which was based on photoionization of 
line-of-sight clumps. However, objections could then arise about justifying the 
complexity of the function found and its overall plausibility. 
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